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ABSTRACT 

The Guitar Nebula is an Ha nebula produced by the interaction of the relativistic wind of a very fast pulsar, 
PSR B2224+65, with the interstellar medium. It consists of a ram-pressure confined bow shock near its head 
and a series of semi-circular bubbles further behind, the two largest of which form the body of the Guitar. We 
present a scenario in which this peculiar morphology is due to instabilities in the back flow from the pulsar 
bow shock. From simulations, these back flows appear similar to jets and their kinetic energy is a large fraction 
of the total energy in the pulsar's relativistic wind. We suggest that, like jets, these flows become unstable 
some distance down-stream, leading to rapid dissipation of the kinetic energy into heat, and the formation of an 
expanding bubble. We show that in this scenario the sizes, velocities, and surface brightnesses of the bubbles 
depend mostly on observables, and that they match roughly what is seen for the Guitar. Similar instabilities 
may account for features seen in other bow shocks. 

Subject headings: pulsars: individual (PSR B2224+65) — ISM: individual (Guitar Nebula) — ISM: bubbles 
— ISM: jets and outflows — instabilities 
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1. INTRODUCTION 

Many pulsars travel at high speed, and the collision between 
their relativistic winds and the interstellar medium leads to the 
formation of bow shocks. These shocks are observed most 
readily at X-ray and optical wavelengths: the shocked rela- 
tivistic wind will emit mostly synchrotron radiation, while the 
shocked interstellar medium will emit - if it is parti ally neutral 
- copious Ha emission (for a review, see iGaensler & Sland 
1200(3) . 

Arguably the most spectacular bow shock is the Guitar Neb- 
ula, made by one of the fastest pulsars known, PSR B2224+65 
(Cordes, Romani, & Lundgren 1993). This Ha nebula has, as 
the name implies, a guitar-like shape, with a bright head, a 
faint ne ck, and a bod y consisting of two larger bubbles (see 
Fig. [TJ. ICordes et al.1 d 19931) suggest this morphology might 
reflect variations in either the pulsar energy injection rate or 
the interstellar medium density. 

In this Letter, we investigate whether instead the peculiar 
morphology could be due to instabilities in the jet-like flow of 
pulsar effluvium away from the bow shock. Fast back flows 
are a natural consequence of bow shocks: the pulsar wind is 
greatly heated at the shock, which, for the usual case where 
cooling is slow, leads to a high pressure that drives a flow in 
the o nly direction available, to the back. From simulations 
(e.g.. [Buccianti ni et alJl2005l) . the flows seem similar to jets, 
being well-collimated and fast, and seem to carry most of the 
pulsar wind energy. Jet-like flows are indeed seen in X-ray 
observations, which also show that only a small fraction of the 
energy is radiated (for a review, Kargaltsev & Pavlov 2008). 

So far, the simulations have not extended far to the back, 
but if the back flow is similar to a jet, one might expect it 
to become unstable further d ownstream. From simulations 
of jets (e.g.. iBodo et al.| [l998). this instability would lead to 
mixing with the ambient medium and rapid dissipation of the 
kinetic energy into heat. The simulations have not followed 
what happens beyond this initial mixing, but it seems plausi- 
ble that the material would expand rapidly and drive a bubble. 
If so, it might initially expand faster than the pulsar motion, 
and gain further energy from the jet. With time, however, it 
will slow down, and once the pulsar has moved sufficiently 



far ahead, the jet will become so long that it becomes un- 
stable before reaching the bubble, and a new bubble will be 
formed. We suggest the body of the guitar is made up of two 
such bubbles, while another one has just started to form near 
the head. 

In §|2j we describe our model in more detail, and in §[3] we 
compare it with the properties of the Guitar Nebula, finding 
qualitative agreement. In §|H we discuss implications as well 
as ways in which our model could be tested. 

2. BLOWING BUBBLES WITH A BOW SHOCK BACK FLOW 

We consider a pulsar that loses energy at a rate E in the 
form of a relativistic wind and moves at velocity v* through 
a medium of density po- The resulting bow shock will have a 
stand-off distance ro given by 
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The bow shock leads to a jet-like back flow carrying kinetic 
energy at a rate fgE, wh ere from simulatio ns the efficiency 
factor /^ is close to unity (Bucciantini 2002). We assume the 
back flow will become unstable some distance £ behind the 
neutron star, rapidly mix with shocked ambient medium and 
dissipate its energy, leading to the formation of a bubble. As- 
suming also that the bubble is fed more energy for some time 
f; n j, and expands adiabatically for a total time f exp , the bub- 
ble radius will be approximately given by the Sedov-Taylor 
solution, 
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where ?7 7 is a dimensionless constant of order unity that de- 
pends on the adiabatic index of the interstellar medium and 
the extent to which energy injection is instantaneous (see be- 
low), and f p takes account of possible variations in density 
between the head of the bow shock and the location of the 
bubble (for our model, by assumption, f p ~ 1). 

Eq.Q]and[2]both depend on the ratio E/p^, suggesting the 
bubble radii can be expressed in terms of bow-shock proper- 
ties and other observables. For this purpose, we rewrite Eq.Q] 
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where i is the inclination, d the distance, ^* = v»sinz'/<i the 
proper motion, and 6$ = firo/d the angular stand-off distance 
(fi i s a function of the inc lination, with /900 = 1, but ft ^sin/; 
see iGaensler et al.l 120021) . We also write f exp = (a-A)//x*, 
where a is the angular separation between the center of the 
bubble and the pulsar and A = isini/d the angular size cor- 
responding to the instability length £, and t m j = (J3- AA)//i*, 
where is the separation between the center of the bubble 
and the next one closer to the pulsar, and AA takes into ac- 
count that two bubbles can have formed at slightly different 
distances behind the pulsar. Note that for the bubble closest to 
the pulsar, f; n j = f exp and one should replace 0— AA with a- A 
below. With this, the angular radius of a bubble, 9\, = Rb/d, is 
given by, 
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where ft, = ?7 7 (/ ^ / f p Jf sin r') 1 / 5 is of order unity. One sees 
that for a bubble far behind the pulsar (i.e., A <C a and 
AA <C 0), there is little room to fiddle: the uncertainties in 
the efficiency factors, geometry, and velocity may amount to 
a factor two, but they enter only to low power. 
With the sizes, the expected Ha photon rates are given by, 
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where f a is the number of Ha photons emitted per neu- 
tral particle before that particle is ionised (~ 0.05 and 0.27 
for case A and B, resp., weakly dependent on velocity; 
IChevalier & Ray mond 197S>), « H o the neutral hydrogen num- 
ber density, v b the expansion rate, and p, b the corresponding 
proper motion, 
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where the coefficient becomes 3/5 for the bubble clos- 
est to the pulsar. Observationally, it is easiest to mea- 
sure surface brightnesses near the limbs of bubbles. For 
measurement length scales 5 <C 0b, one predicts s ay n mb = 

n a . b {\/2n5e h )^25j¥h^^X 12 - 
We now discuss our assumptions and simplifications. Our 

three main premises are that the back flow is jet-like; that 
it becomes unstable; and that the instability drives a bubble 
that is fed further energy for some time. The first two are 
supported by simulations: bow shocks appear to give jet-like 
back flows with rj et ~ 4ro (Buccia ntiniet alj 2005), and jets 
do seem to become unstable (with much of the recent work 
focussing on how to prevent th is from happening too quickly; 
for a review, e.g.. lHardeell2004l) . Typically, perturbations ap- 
pear to grow on length scales £ g about ten times the jet ra- 
dius rjet. For instance, f or relativistic jets in a ten times denser 
medium, lHardee et all d!998l) found 4 < L/r iet < 15. For a 
bow shock, this implies angular growth length scales A g in the 
range 16 < (X g /6o)(Ji/smi) < 60. Of course, instability will 
only occur after a few growth times, so the angular distance A 
should be correspondingly larger. We will see in §[3]that this is 
consistent with the Guitar Nebula. It also validates an implicit 
assumption we made, that the bubbles do not overtake the pul- 
sar, i.e., that a > 8 b at all times. From Eq. [4] the minimum 



value of a- 0b occurs at a = A + 6»o(3/ b /5) 5/2 (47rc/v») 1/2 ; this 
becomes negative only for v* < 50/^(50#o/A) 2 kms -1 , much 
smaller than the velocity of PSR B2224+65. 

What is not clear yet, however, is whether jet instabili- 
ties could lead to bubbles, and, if so, whether our simpli- 
fied description is justified. In particular, in using the Sedov- 
Taylor solution, we assume energy is injected (nearly) in- 
stantaneously at a single point in a homogeneous medium, 
and that the expansion is adiabatic. Of these assump- 
tions, the last is reasonable: the cooling time, f coo i w 
^OOOO yrtvb/lOOkms-'^po/Kr^gcm- 3 )- 1 (IK00 & McKed 
119921) . is much longer than the ~ 300 yr it takes PSR 
B2224+65 to cross the nebula (here, we scaled to the low- 
est velocities and highest densities appropriate for the Guitar 
Nebula). The others are less realistic: energy will be injected 
some time in a larger, not necessarily spherical volume em- 
bedded in a medium which, close to the axis, has been through 
the bow shock. Injection over some time should lead to a 
bubble that initiall y expands somewhat more slowly. Indeed, 
Dokuchaev (2002) found that for continuous energy injection, 
Eq. [2] can be used, but with a somewhat smaller value of 7? 7 
(e.g., 775/3 = 0.929 for f exp = t m ^ instead of 775/3 = 1.152 for 
finj <C f eX p)- This may be counteracted, however, by the ini- 
tial expansion being in pre-shocked, less dense medium. At 
later times, our estimates should depend less on these initial 
conditions, but rather on the extent to which bubbles can be 
treated in isolation, when in the Guitar Nebula they appear to 
have merged (Fig. [TJ. Overall, we conlude that our heuristic 
model will only be good at the factor 2 level. 

3. RECONSTRUCTING THE GUITAR 

To see how well our model applies to the Guitar Nebula, 
we retrieved a deep Ha im age taken on 6 Dec. 2000 for the 
Champlane survey dZhao et al.l 2005), and measured proper- 
ties of what seemed the three most obvious bubbles (see Fig.Q] 
and Table [TJ: the two forming the body of the guitar, and one 
just behind the pulsar (hereafter, bottom, middle, and head). 
For all bubbles, the estimates of the angular separation from 
the pulsar, a, and of the angular radius, 9 b are quite reliable, 
but for the middle bubble, the separation to the next closest 
bubble, 0, is relatively poorly defined, since there may be an 
additional bubble in the neck. 

F or our e stimates, we also need the stand-off distance 
6» . Chatt erjee & Co rdes (2004) find that the shape of the 
bow shock - as seen in Ha images taken with the Hubble 
Space Telesco pe (HST ) - is reproduced well with the ana- 
lytic model of I Wilkin! (119961) . For data sets taken in 1994 
and 2001, they infer inclinations i ~ 90° and stand-off radii 
a = Of 12 ±0"04 and 0."15 ±0"04, respectively. Since Ha is 
emitted outside the actual stand-off distances, with 6 a ~ 1.5#o 
dBucciantinil |2002|) . one infers 9 Q ~ 0."09, which we will 
use below. We wi ll also use the obser ved proper motion 
of ft* = 0"182yr _1 (lHarrison et al.lll993l) and scale the pul- 
sar speed to v* = 1500 km s" 1 (where we used i ~ 90° and a 
distance d = 1.8kpc, as implied by the dispersion measure 
of 35.30pccm~ 3 and the NE2001 electron density model of 
ICordes & LazToll2002l) . 

3.1. Sizes 

In matching our model to the measurements, we first note 
that the existence of the head bubble implies A < 9". Thus, 
for the bottom bubble, A <C a and AA < A < 0. For this 
case, Eq.Hsimplifies to 6> b =/ b ((47rc/v»)6» 2 a 2 /3) 1 / 5 = 18"3/ b ; 




TABLE 1 

Bubbles in the Guitar Nebula 



20 40 60 

z (arcsec) 

FIG. 1. — Reconstructing the Guitar. (Top) A possible sequence of events 
that could have led to the guitar nebula, using instabilities in the jet-like pulsar 
back flow. The instabilities are set to occur at the positions marked with 
a star, roughly 8" behind the pulsar. They lead to a rapid dissipation of 
energy and spherical expansion, and are fed for a while, until an instability 
occurs closer to the pulsar. To match roughly the image of the guitar, we 
have to assume that the instabilities occurred rather closely together in the 
region marked with < signs, so that no distinct bubbles formed. (Bottom) 
Continuum-subtracted Ha image from the Champlane survey, with suggested 
locations of bubbles overlaid. 

to match the observed radius of 16" thus requires fi, ~ 0.9, 
close to unity as expected. 

Using this value of fi, for the head bubble, we find we re- 
quire A ~ 7."5 to match the observed small size of 2". This 
implies X/9(, ~ 80, in line with expectations (§|2). It also im- 
plies the bubble formed only recently, about (a- A)//i» ~ 8yr 
before the image was taken. We return to this below. 

For the middle bubble, we find that to match its observed 
size requires j3- A A ~ 5 "5. This is much smaller than for the 
bottom bubble, since to produce this relatively small bubble 
requires much less energy. It raises the question, though, what 
happened to the energy dissipated later, outside the middle 
bubble. One possibility is that more energy was injected in 
the middle bubble, but that it grew, it merged and equilibrated 
with the larger bottom bubble. If so, our above estimate of 
/b would be too large. Clearly, we have reached the limits of 
applicability of our simplistic picture of individual, spherical 
bubbles. 

3.2. Proper motions and brightnesses 

Our model predicts expansion rates (see Table [TJ. For the 
head bubble, the predicted rat e is fast, ^O^lSyr " 1 . By com- 
paring HST images, Chatterie e & Cordesl (120041) indeed find 
that the head bubble expanded between 1994 and 2001, es- 
pecially to the back, at a rate comparable if slightly slower 
than that predicted, of ~0."10yr _1 . Interestingly, the bubble 
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NOTE. — For each bubble, the first row lists measured parameters: 
o, the angular separation between the centre of the bubbles and the 
pulsar; 9, the angular radius; j3, the separation to the next bubble (un- 
defined for the head bubble); fi, the expansion rate; and s a ji m b, the 
surface brightness near the limb (relative to that at the position of the 
pulsar, which, from the instrumental sensitivities, has a photon rate of 
1.0 X lQT* s cnT 1 arcsec -2 , with an uncertainty of about 20%). The 
second row lists model assumptions and predictions, including A, the 
angular wavelength of instability. Empty entries indicate that observed 
values were used. Limb surface brightnesses for the middle and bot- 
tom bubble were calculated relative to the observed brightness for the 
head bubble, using that S a Aimb °c p9 1 ' 2 (see §f2). 

also became brighter, consistent with the idea that it formed 
only recently. This evolution is confirmed by inspection of 
unpublished HST data taken in 2006. Furthermo re, the head 
bubble is dimmer in the 1992 discovery image of Corde s et alj 
dl993l) than it is in Fig. [1] or in the 1995 image shown by 
IChatteriee & Cordesl d2002l) . 

For the middle and bottom bubbles, the predicted expansion 
rates are slower, ^0."02yr _1 . This is difficult to detect from 
the ground. It may be detectable over the 12 years spanned by 
the HST images, but given the low signal-to-noise ratio, this 
will require detailed modelling, which we have not attempted. 

The lower proper motions for the middle and bottom bub- 
bles also imply predicted limb surface brightnesses about 3 
times fainter than for the head bubble. This is roughly consis- 
tent with the observed ratio of 4 (Table[T|i. 

4. RAMIFICATIONS 

We found that we could roughly reproduce the Guitar Neb- 
ula assuming the jet-like back flow from the pulsar bow shock 
becomes unstable and dissipates rapidly, causing expanding 
bubbles. If this were to happen generally, one might expect 
other sources with jets or bow shocks to show Guitar-like bub- 
bles, yet none appear to be known. For jet sources, this may 
not be surprising: many jets are denser than the medium they 
move through, and hence more stable, and disruptions that do 
occur may be difficult to distinguish from, e.g., changes in jet 
orientation. 

For other bow shocks, the absence of bubbles may partly be 
a selection effect: most have much larger stand-off radii than 
the Guitar, and hence any bubbles would be at correspond- 
ingly larger distances, where they might be missed, especially 
as they would be fainter than the bow shock (or even invisible 
if the expansion velocity became too low or if radiative ef- 
fects became important; both perhaps relevant especially for 
stellar wind bow shocks). The one possible exception is PSR 
B0740-28, which has a Ha bow shock with a relatively small 
stand- off radius of flp = 1 "0 as well as "shoulders" further be- 
hind d Jones et al.ll2002l) . If related to an instability, one infers 
15 < A/#o ;$ 60, of the same order as we see for the Guitar. 
It would be interesting to obtain deeper images further behind 
the bow shock. 

In some pulsar bow shocks, the shocked pulsar wind is ob- 



served directly, by its syn chrotron emission (for an overview, 
Kargaltsev & Pavlov 2008). For many, including the Gui- 
tar dHui & Beckerll2007bl) , emission is seen only close to the 
pulsar, likely at the pulsar wind termination shock. Some, 
however, have much longer tails. The longest belongs to 
the "Mouse," associated with PSR J1747-2958. This neb- 
ula, with #o ~0."75, shows a bulbous structure ~ 1.'5 behind 
the pulsar (the Mouse's body ), but also a smo oth, straight tail 
of 12', without a clear end dGaensler et al.l 2004). Scaling 
with the stand-off distances, one might identify the Mouse's 
body with the equivalent of the Guitar's head bubble. The 
long tail has a size equivalent to the bottom bubble, but, 
apart from changes in pol arisation, shows little structure 
dYusef-Zadeh & Balrvlll987l) . This would seem inconsistent 
with any bubbles being formed, and thus is puzzling in the 
context of our model. 

For two other pulsar bow shocks with long tails, the ob- 
servations match expectations better. For PSR J1509-5850, 
with 6>o ~ 0"5, the X-ray tail extends for >5.'6 and shows 
clear structure, with a change in brightness at 1.'3, a kink 
at 3', and a bright radio spot coincident with its end point 
dHui & Becker! l2007al; iKarealtsev et al.l 120081) . Comparing 
with the large bubbles in the Guitar, the typical length scale of 
~ 1 '.5 for the knots and kinks is about a factor 3 larger, roughly 
consistent with the ratio of the stand-off distances. For PSR 
B1929+10, with 9 ^ 2"3, the tail extends up to 10' and again 
shows substantial structure, with brighteni ngs at ~2' and ~5', 
the latter coincident w ith a radio feature dBecker et al.ll2006t 
Misan ovic et al.l2007l) . Again scaling with the stand-off radii, 
the 5' feature could be similar to the head bubble in the Guitar. 

Overall, we conclude that our model of instabilities in a 
bow shock back flow roughly reproduces observations of the 



Guitar Nebula, without the need to appeal to variations in the 
density of the ambient medium, nor to energy sources beyond 
what is expected to be carried by the back flow. It also seems 
consistent with what is seen in other pulsar bow shocks. The 
model could be tested further both with observations and sim- 
ulations. Observationally, one test would be to measure the 
expansion velocities in the Guitar bubbles, either by deter- 
mining proper motions, or by spectroscopy (from the broad 
component of the Ha p rofile, as done for non-radiative shocks 
in supernova remnants; Raymond 1991). Given the observed 
Ha surface brightness, this would allow one to estimate the 
ambient density, which should be similar to that at the loca- 
tion of the bow shock in our model, but substantially lower 
if the bubbles reflect dens ity variations dCordesetal .11 19931; 
IChatteriee & Cordesll20"04T) . 

Simulations of bow shocks that extend to larger scales 
might show whether instabilities in fact lead to bubbles or 
rather to more continuous structure, or whether perhaps the 
process is sufficiently stochastic that both can occur (possibly 
leading to a shape like the Guitar's neck). If bubbles form, the 
simulations might also shed light on details of the morphol- 
ogy, such as the closed appearance at the back of the head and 
bottom bubbles. 
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